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We present a novel system to calibrate millimeter-wave polarimeters for CMB polarization measurements. 
This technique is an extension of the conventional metal mirror rotation approach, however it employs 
cryogenically-cooled blackbody absorbers. The primary advantage of this system is that it can generate 
a slightly polarized signal 100 mK) in the laboratory; this is at a similar level to that measured by ground- 
based CMB polarization experiments observing a ^ 10 K sky. It is important to reproduce the observing 
condition in the laboratry for reliable characterization of polarimeters before deployment. In this paper, we 
present the design and principle of the system, and demonstrate its use with a coherent-type polarimeter 
used for an actual CMB polarization experiment. This technique can also be applied to incoherent-type 
polarimeters and it is very promising for the next-generation CMB polarization experiments. 
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I. INTRODUCTION 

The faint pattern of cosmic microwave back- 
ground (CMB) polarization promises to provide new and 
interesting information on the universe. The primary 
goal of CMB polarization studies in the next decade is 
to detect the degree-scale B-modes (curl components) in- 
duced by primordial gravitational waves. The existence 
of the primordial gravitational waves is a generic predic- 
tion of inflation. Therefore the detection of the i?-modes 
is a "smoking-gun" signature of inflation^. Using an ar- 
ray with ~ 1,000 of polarization detectors (hereafter po- 
larimeters) is essential to discover the _B-modes^. For the 
next-generation experiments with such a large polarime- 
ter array, establishing the performance of the detectors 
in the laboratory is essential to prepare for observations 
in the field. 

Our system aims to calibrate three important prop- 
erties; responsivity, orientation of detector polarization 
angle and spurious polarization signal in the instrument. 
The calibration of these parameters requires a well- 
characterized artificial polarization signal. The combi- 
nation of an unpolarized cold load (blackbody absorber) 
and a metal mirror at room temperature is one of the 
main approaches^r— to obtain such a signal. The temper- 
ature difference between the cold load and the mirror cre- 
ates a well-characterized polarization signal, ~ 100 mK 
(Details are discussed in Sec. Ill A[) . In the laboratory, 
the blackbody absorber tends to be cooled with liquid 
nitrogen (77 K). However, the large difference between 
this temperature and that of the sky at the observing 
site (^ 10 K, for example the Atacama desert of north- 
ern Chile at an altitude of ^ 5000 m) makes it difficult 
to characterize the polarimeters in the laboratory. This 



is particularly true for detectors with a narrow range of 
linear response to input load temperature. 

To address these issues, we propose an advanced ap- 
proach for polarimeter calibration. Our method is an 
extension of the metal mirror rotation approacb^r^, how- 
ever it employs cryogenically-cooled blackbody absorbers 
instead of the absorbers cooled with liquid nitrogen. This 
system provides in the laboratory a load condition similar 
to the actual observation. 



II. PRINCIPLE OF POLARIZATION CALIBRATOR 

A. Principle of polarization signal generation 

The principle for generating the polarization signal is 
the same as the conventional method'^"^ . Figure [T] shows 
a schematic view of the developed calibrator system. The 
blackbody emitters are cooled with a cryocooler as the 
cold load with temperature of Tload- Unpolarized black- 
body radiation is emitted from the cold load; it reflects 
off a metal mirror surface, which induces a linearly po- 
larized component because of the finite emissivity of the 
mirror. The magnitude of the polarization signal (P) can 
be calculated using the following formula. 



P = 2y/l6Tn^peo tan {(3) {T,,- 



Tload) , (1) 



where v is the frequency of observed radiation, Tmirror 
and p are the temperature and resistivity of the mirror, 
/? is the refiection angle of the radiation, and eq is the 
vacuum permittivity. For example, we obtain ^110 mK 
polarization signal when we use an aluminum mirror (p ^ 
2.2 pn ■ cm) with Tmhror = 250 K, Tioad = 10 K and /3 = 
45°. We can control the signal magnitude by changing p 
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FIG. 1. Schematic view of tlie developed calibration system. Unpolarized radiation is emitted from the blackbody absorbers 
cooled with a cryocooler. The radiation reflects off a metal mirror surface placed at the center, and enters a tested polarimeter. 
The incident radiation is polarized because of the finite emissivity of the mirror. The blackbody absorbers are cylindrically 
placed as shown to the right. The cryocooler in the system has two cooling stages, and a part of absorbers are cooled to be 
different (higher) temperature using another cooling stage. The mirror can rotate in the vacuum, and the loads at two different 
temperatures are alternately presented to the polarimeter. 



as well as changing the temperature difference, Tmiirc 

Tload- 



B. Principle of Polarimeter Calibration 

The linear polarization is characterized by two Stokes 
parameters Q and U^. Since the value of Stokes pa- 
rameters depends on the bases of the selected coordinate 
system, the ability to rotate the axis of the polarization 
signal is essential to calibrate the polarimeter responses 
unless we calibrate the orientation of the polarization an- 
gle for the Stokes parameters. 

Moreover, the baseline of the polarimeter response fluc- 
tuates due to a 1// noise. To measure the response for 
the polarization correctly, we have to modulate the polar- 
ization signal by rotating the mirror; the parallactic angle 
of the signal axis with respect to the polarimeter varies 
with the rotation angle of the mirror {6). For example, 
assuming the axis for the Q response of the polarime- 
ter aligns with the orientation of 6* = 0, the Q response 
varies with P cos (26'), and the U response varies with 
Psin(26l) (=Pcos(2(6l- f))). We can characterize the 
responsivity; signal height with respect to the input po- 
larization signal, and the polarization angle; orientation 
angle of the polarization axis from the amplitude and 
the phase of the sinusoidal response of the polarimeter, 
respectively. 

We are also interested in the spurious polarization of 
the polarimeter with respect to the total power which is 
usually dominated by unpolarized temperature. This is 
called I to Q (or U) leakage where the total power is de- 
scribed with Stokes / parameter. The actual polarization 



response for Q is Pcos(20)-f e/, where e is the leakage pa- 
rameter Since / is dominated by Tioad in most cases, we 
can measure the leakage parameter for the polarimeter 
by changing the Tioad^. With the two different tempera- 
ture loads of 10 K and 30 K, we observe a 1 K baseline 
difference in the polarization response for e — 0.05. 

The polarimeter can also measure the unpolarized tem- 
perature response for / as well as for Q and U . We cal- 
ibrate the response for / by using two different Tioad s 
(the Y- factor measurement). This allows us to measure 
receiver temperature (Tj-ec) which is the effective temper- 
ature offset from the intrinsic noise. 



C. Expected Polarimeter Response 

In general, the polarization angle does not align with 
the design value perfectly. The responses of the polarime- 
ter are defined as described in the following formula in 
the unit of mV, 



Vq = Rq [P cos (2(0^ 

+eQ(rioad + P) ], (2) 
Vu = Ru [Psin(2(0-0c/)) 

+ec/ (Tioad + P) ] , (3) 

Vl = Rl [ Tioad + P + Tlec ] , (4) 

where Rk {k = Q, U) and Rj are the responsivities 
(mV/K), (f>k is the polarization angle and Sk is the leak- 
age parameter for each Stokes parameter response. 

Figure [2] shows the polarimeter responses as a func- 
tion of the rotation angle of the mirror. Here, we 
assume a configuration with two different temperature 
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.2. i he expected response oi total power and polarization 
signal plotted as a function of the rotation angle (6). The two 
different temperature pyramid arrays are alternately viewed 
by the rotating metal mirror. 



FIG. 3. CR-112 blackbody emitter array as the cold load. 



loads Tt,^ 



10 K and T",, 



rameter e = 0.05. The T^^^^ 



= 30 K and the leakage pa- 
load covers the 10° - 300° 



region and other regions are covered by the T^^^ load. 
The boundary regions between the two loads are masked 
in the analysis because it is complicated to estimate the 
response due to the beam profile. 

We can extract all the calibration parameters from the 
simultaneous fit for each polarimeter response with re- 
spect to the mirror rotation angle (9) and the load tem- 
peratures (Tioad)- In the case Tioad is constant, the po- 
larimeter responses Vk are simple sinusoidal shape with 
constant offsets. Furthermore, in the case CfcATioad ^ P 
(where ATioad = T^ad ~ "^load)' leakage parameters 
£k can be extracted from the baseline difference A 14 be- 
tween ^ and TL^d > i-e- £fc = AH/(i?fcArioad). 
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FIG. 4. Measured power reflection coefficient of CR-f 12 and 
CV-3. 



III. SYSTEM COMPONENTS 



A. Cryostat and Cryocooler 

A cylindrical cryostat with a dimension of 540 mm di- 
ameter and 500 mm height houses all of the system com- 
ponents as shown in Fig. [TJ The system is cooled with 
the 2-stage Gifford-McMahon cryocooler^ Its cooling ca- 
pacity in the cryostat is consistent with the specifica- 
tion (Appendix 1^ . Thermal load for the 1st (2nd) stage 
in our system is about 18 W (7 W), and we usually op- 
erate at 30 K (10 K) temperature on the top of the 1st 
(2nd) stage. 



B. Cold Load - Blackbody Emitter Array 

To create the low temperature blackbody radiation, 
low reflectivity and high thermal conductivity are re- 
quired for the blackbody emitters. Furthermore, the 
emitters need to have viscosity (or similar thermal ex- 
pansion coefficients with the base plate material) to with- 



stand the thermal cycling. To fulfill the requirements, we 
employ Eccosorb CR-112^ an iron-loaded epoxy as the 
material for the blackbody emitter. To minimize the sur- 
face reflection, the CR-112 is casted with pyramid-shape 
grooves on the front surface. We follow the manufactur- 
ing processes of the ARCADE experiment^. The pyramid 
has a square base of 400 mm^ and height of 58 mm. To 
achieve higher thermal conductivity and uniform surface 
temperature, each pyramid has an aluminum core and a 
copper wire epoxied onto the end of the aluminum core, 
running almost to the tip. We use two arrays of the pyra- 
mids as the cold load (Fig. [3]). We use 420 pyramids in 
total to surround the inside of the cryostat wall in the 
system. 

We confirmed the low refiectivity of the pyramid ar- 
ray by measurements with commercial millimeter-wave 
equipmentiiS. The measured reflectivity for the pyramid 
array is shown in Fig. 01 It is about —30 dB around the 
90 GHz frequency region and is comparable to commer- 
cial absorber material such as CV-3*. 
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C. Metal Mirror 

We adopt /3 — 45° for the metal mirror in our system. 
The mirror is hung with a stainless steel shaft that is 
rotated with a DC motor from the outside of the cryo- 
stat via a feedthrough^^ The maximum rotation speed is 
~12 rpm, which corresponds to much higher frequency 
of the 1// noise of the polarimeters for the CMB exper- 
iment (several mHz ~' 100 niHz). 

We use three metal mirrors, aluminum, steel and stain- 
less steel, to vary the polarization signal intensity. It is 
useful to evaluate the possible inherent polarization^ii^. 
We obtained zero-consistent inherent polarization in our 
system (Sec HVD]) . 

Currently, precision of the resistivity (p) is the main 
uncertainty for the calculation of the polarization signal 
(P). The variation of the catalog specification is 20%, 
corresponding to a 10% uncertainty for the polarization 
signal calculation. Such uncertainty can be reduced with 
a resistivity measurement in a separate setup. 



D. Load and Metal Mirror Temperature 
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FIG. 5. Loads and mirror temperatures during the cooling 
process. 



Figure [S] shows the temperature of the loads and the 
aluminum mirror. It takes ~10 (~ 30) hours to stabi- 
lize the temperature of the cold load (mirror) Because 
of the difference of the emissivities for each metal, equi- 
librium temperatures depend on the mirror material as 
summarized in Table HI The measured temperatures are 
consistent with the thermal calculation in Appendix |B] 

Uniformity of the load surface temperature is con- 
firmed to be better than 1 K with DT-670 commer- 
cial silicon diode thermometers^^ The temperature non- 
uniformity causes only a 0.4% effect for the polarization 
signal while the effect for the unpolarized response is 
about 10%. The temperature uniformity can also be con- 
strained from the flatness of the Vj response as a function 
of the rotation angle effectively, which is also less than 
1 K fSec HVB]) . 



IV. EVALUATION OF SYSTEM PERFORMANCE 

A. Polarimeter for System Evaluation 

We evaluate the system performance using the QUIET 
polarimeter—. QUIET has an array of pseudo-correlation 
polarimeters, each of them employs hybrid couplers to 
correlate the two orthogonal polarizations selected by an 
orthomode transducer in order to directly and simultane- 
ously measure the Stokes Q, U and I parameters of the 
incoming signal. A corrugated horn is attached on top of 
each polarimeter, and its beam width is measured to be 
^ 10° (FWIIM)^'^^ The beam coverage of this system is 
^ 20° , and the spill-over power is expected to be less than 
1%. The responsivity of this polarimeter is ~ 0.15 mV/K 
including the pre-amplifier gain factor of ~ 100^^ which 
was measured with a different setup using the rotation 
plate technique under 77 K. In case of QUIET polarime- 
ters, reflections around the septum polarizer is the main 
reason for the I to Q [U) leakage. A fraction of the 
power input on one port of the correlation module was 
reflected to the septum polarizer and a fraction of the 
reflected power re-entered the other port. Typically, the 
leakage of a QUIET polarimeter is £ 0.005. To evalu- 
ate the system performance for the leakage parameter (e) 
easily, we insert an extra waveguide between the corru- 
gated horn and the septum polarizer. This increases the 
reflection, and the leakage size is increased to e ^ 0.05. 



B. Total power response 

The response for the total power as a function of the 
rotation angle of the mirror is shown in the left panel 
of Fig. [6] The total power response shows two levels; 
the upper (lower) level corresponds to 10 K (30 K) load 
temperature. The total power responsivity, i?/, is calcu- 
lated from the timestream with Eq.|3]to be 0.15 (mV/K), 
which is consistent with the previous measurement dis- 
cussed in Sec. IIV Al 

The uncertainty for the total power response is domi- 
nated by the uniformity of the load temperature, which 
is evaluated from the timestream of the total power re- 
sponse. The RMS of the total power response for the 
^load lo^'i is 0.9 K, which is consistent with the RMS 
on the temperature monitored using thermometers. This 
variation includes the polarimeter noise which dominates 
the response fluctuation in this test. Thus, this number 
still provides an upper limit on the non-uniformity of 
the load temperature. So far, the precision for the total 
power responsivity measurement is 7% (1 K uncertainty 
both for Tj^jjj and T^^^, which is already comparable 
with the precision obtained using astronomical sources. 
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TABLE I. Resistivity (p), load temperatures and mirror temperature for each mirror material. P is the calculated polarization 
amplitude which corresponds to each T^^^. In case of the stainless steel mirror, we use an aluminum shaft, which has better 
heat conductivity than that of a stainless steel shaft. Therefore Tmirror is higher than others. 



Mirror Material 
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FIG. 6. Response of the QUIET polarimeter to the rotating aluminum block. The left panel shows the timestream for a total 
power measurement and the right panel shows the timestream for a polarization measurement. A red dashed line in the right 
panel is a fitted sinusoidal curve. 



C. Polarization response 

The response for the polarization as a function of rota- 
tion angle of the aluminum mirror is shown in the right 
panel of Fig. [51 The response shows a clear sinusoidal 
curve as a function of the mirror rotation angle. This is 
the first demonstration "in the laboratory" for the mod- 
ulation of the polarization signal (110 mK) with a load 
temperature which is as low as those at the observation 
site. 

From the amplitude of the curve, we obtain Rjj = 
0.12 ± O.Ol(siai) ± Qm[syst) mV/K which is consistent 
with the i?/. Here the stat term is the statistical un- 
certainty due to the polarimeter noise and the syst term 
is the systematic uncertainty derived from the catalog 
uncertainty of the resistivity (10%, Sec. IIII C[) and the 
load temperature (0.4%). This is the uncertainty for the 
absolute scale, not for the relative uncertainty among de- 
tectors. The relative uncertainty is much more important 
for CMB polarimeters. 

For the polarization angle, (pu, we extract the angle 
with a 0.8° error from the fit, where the error is dom- 
inated by the polarimeter noise. The system has not 
been limited by the precision of the relative angle mea- 
surement. However, we have not implemented the en- 
coder for the mirror rotation angle yet, which limits the 
absolute angle measurement. This encoder is needed to 
achieve the necessary precision of 0.5° for future i?-modes 
measurementaii. 

For the I to Q (or U) leakage, we can determine the 
leakage parameter with a precision of 0.005 from Fig. [51 



Again this measurement is limited by the polarimeter 
noise. There is no systematic limitation for the leakage 
measurement unless the relative uncertainty between i?/ 
and Rq^u becomes significantly large. 

Just a few tens of seconds of measurements using our 
system can provide comparable precision with astronomi- 
cal source calibration at the observation site, for example 
a 20-minute observation for the Crab nebula (Tau A)^-. 
Our system therefore provides better precision than as- 
tronomical sources in a given amount of time. 



D. Evaluation of Inherent polarization signal in the system 

Figure[7]shows the amplitude of the polarized signal for 
three different mirror materials. The results show that 
the amplitudes are proportional to the square root of the 
resistivity as expected, and the system has no inherent 
polarized signal within the precision of the measurement. 
This is the evidence that there is no bias for the polariza- 
tion rcsponsivity measurement by using the single mirror 
setup. 



V. DISCUSSION AND SUMMARY 

We have developed a polarization calibrator for CMB 
polarization receivers. The technique is a simple exten- 
sion of the conventional rotation metal plate approach, 
however, it employs ~10 K cold loads which correspond 
to the sky temperature at the CMB observation site 
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FIG. 7. The response of the polarized signal is plotted with 
respect to the square root of the resistivity of the mirror mate- 
rial. Here three different materials (Aluminum (•), Steel (■) 
and Stainless steel (*)) are used. 



on the ground (the Atacama desert of northern Chile 
at an altitude of ~ 5000 m). The calibrator can pro- 
vide a well-characterized polarization signal (^ 100 niK) 
which can be modulated by rotation of a mirror. Us- 
ing a QUIET polarimeter, we demonstrated simultane- 
ous measurements of the responsivity, orientation of po- 
larization angle, and spurious polarization signal in the 
instrument with sufficient precision. These three param- 
eters are keys to the success of B-mode observations as 
well as sensitivity (noise equivalent temperature: fxK^/s) 
of the polarimeter. Our system also provides a way to 
measure noise properties under a comparable load tem- 
perature condition to the observing site. The estimation 
based on our system is thus more "reliable" than that ob- 
tained using a 77 K load cooled with liquid nitrogen. Our 
system can easily be extended for large detector arrays 
by scaling up the size of the system. 

The precision of the laboratory-based calibration is 
comparable with the astronomical source calibrations 
at the observation site. This system can establish the 
polarimeter performance before deployment, which is 
important for experiments with large detector arrays. 
Moreover, the precision of calibrations using astronom- 
ical sources is limited by the beam profile uncertainty. 
Our system is free from this effect because of the wide 
coverage-angle of the loads. We have already eliminated 
it to less than 1% for the QUIET polarimeter. 

We have confirmed that our system works as designed. 
It will be very useful for the next-generation CMB polar- 
ization experiments. In this paper, we demonstrated the 
performance of the system using a polarimeter based on 
coherent amplifiers. By using a different configuration, 
e.g. a setup in which the refiected emission on the mirror 
goes outside the cryostat, it is possible to use the system 
as an external polarized source for a polarimeter with 
different technology such as a bolometer. With a larger- 
capacity cryocooler, the system can also be extended to 
3 K load temperatures. Such a system will be a key tool 



in the development of polarimeters for the future satellite 
experiments. 
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Appendix A: Evaluation of Cooling Capacity of Cryocooler 





20 






III 


1 1 1 1 1 1 1 1 




18 










a 










3 


16 










^-^ 












re 

0) 


14 




(13.0 W, 0.0 W) 




(13.0 W, 35.4 W) 


mp 


12 










0) 
♦J 


10 




(7.0 W, O.OW)^^^- -■- 







0) 


8 






Operation 




O) 








re 










^-^ 
<n 


6 










•a 












c 
o 


4 




(0.0 W, 0.0 W) (0.0 W, 1 


8.0 W) 


(0.0 W, 35.4 W) 


o 


2 

























15 20 25 30 35 40 45 
First stage temperature [K] 

FIG. 8. Measured Load map of the cryocooler used in the 
system, SHI-RDK 408S, which is consistent with the sepcifi- 
cation. 



Figure [8] shows the measured load map for the 10 K 
Gifford-McMahon cryocooler (Sffl-RDK 408S) used in 
the calibrator system. The first (second) stage of the re- 
frigerator reaches temperature of 22 K (4.6 K) under no 
external thermal load. The calibration system is usually 
operated at 30 K (10 K) on the first (second) stage. The 
thermal load on each stage during the operation is esti- 
mated to be 18 W (7 W) from the load map measurement. 
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FIG. 9. Schematic diagram of a cryostat for lieat transfer 
calculation on the mirror. 



Appendix B: Heat transfer on the metal mirror and mirror 
temperature in the system 

Figure [S] shows a schematic diagram of the heat trans- 
fer on the metal mirror. The mirror has a surface area of 
0.04 m^ and is supported by a stainless steel shaft (the 
surface area divided by the length is 0.09 m). The mirror 
temperature instationary state is estimated as the tem- 
perature where the following heat transfer balances. 

1- Qrad : Radiant heat inflow from the cold loads and 
outflow. 

2. Ocon : Heat conduction from a DC motor through 
the shaft. 

Each heat transfer is calculated as follows; 



Qrad — f ^11 



A 



V mirror 



L 



300K 



k{T)dT 



^load) £mirror (Bl) 

(B2) 



where cr is Stefan's constant (5.67 x 10"^Wm"^K-^), 
^mirror IS the surface area of the mirror, Tmirror and Tioad 
are the physical temperature of the mirror and the cold 
load, emiiTor is the emissivity of the mirror, A/L \s the 
ratio of the surface area to the length of the stainless 
steel shaft, and k{T) is the temperature-dependent ther- 
mal conductivity of the shaft, respectively. In case we 



use an aluminum mirror (^mirror 0.03 ~ 0.04), the heat 
transfers are balanced at Tminor of 240 K ~ 260 K. This 
is consistent with the measured temperature as stated in 
SeciniDl 
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